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I. INTRODUCTION 

The static correlations between charged particles in classical plasmas at equilibrium 
have attracted some theoretical attention [1] . For a plasma with a boundary, the surface 
charge correlations are especially interesting, because they are universal, in that sense that 
they do not depend on the detail of the microscopic constitution of the plasma, for length 
scales large compared to the microscopic ones [2] (although these correlations do depend 
on the geometry of the boundary). 

Experimentally [3,4], classical plasmas in thermal equilibrium, with a boundary, have 
been obtained by confining particles of one sign in a Penning or Paul trap. These particles 
may be electrons or ions. In a Penning trap, a magnetic field along the z axis confines 
the particules radially while an electric field due to suitable electrodes provides the axial 
confinement. In a Paul trap, the confinement is provided by an electrostatic potential 
and a radiofrequency field. In both cases, one obtains a system which behaves like a one- 
component plasma: a system of particles of one sign immersed in a neutralizing uniform 
background (in the experiments, the confining fields play the role of the background). The 
particles form a blob of spheroidal shape. In the case of a Penning trap, the blob performs a 
rigid rotation around the z axis, and it is in the corotating frame that it behaves like a static 
one-component plasma in a magnetic field; however, if the rotation frequency is just half 
the cyclotron frequency (the Brillouin regime), the system behaves like a one-component 
plasma without a magnetic field. 

Electrostatic collective modes of such spheroidal plasmas have been experimentally 
observed [5-8] and theoretically discussed [9-11]. 

These experimental findings are an incentive for applying the general theory [2] of 
surface charge correlations to spheroidal plasmas; this will be done in Section II. Further- 
more, since there is some knowledge about the electrostatic collective modes, it is desirable 
to check in a variety of cases that the static charge correlations can be correctly obtained 
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(on large length scales) as the sum of the contributions from all the thermally excited col- 
lective modes (it was found some time ago that this approach does work in the very simple 
special case of a plane boundary without a magnetic field [12]); this will be done in Sec- 
tion III. This method for computing the static correlations also gives the time-dependent 
correlations. 

II. STATIC CORRELATIONS IN A SPHEROIDAL PLASMA 

It should be emphasized that the calculations of this Section are based on macroscopic 
electrostatics and give results which are valid only on macroscopic length scales. There- 
fore, these calculations cannot account for the shell or crystal structures which have been 
observed [13] for plasmas in a Penning trap and reproduced by computer simulations [14]. 
The present calculations apply to correlations smoothed on microscopic oscillations. Also, 
the surface charge density a which will be considered here should be understood as a mi- 
croscopic volume charge density integrated along the normal to the surface and smoothed 
in directions parallel to the surface. 

As explained in the Introduction, a plasma confined in a Penning trap behaves, in 
a rotating frame, like a static spheroidal one-component plasma submitted to a uniform 
magnetic field along its axis. We assume that the Debye length is much smaller than the 
plasma dimensions: this condition ensures that the spheroid has a well-defined surface and 
that the spheroid size is macroscopic. 

It is well known that in classical statistical mechanics a magnetic field has no effect 
on the static quantities. Thus we shall compute the static (equal-time) correlations of a 
spheroidal one-component plasma in the absence of a magnetic field. The results will be 
also applicable to the case with a magnetic field. 

We use the method of Ref. [2] . Surface charge correlations can be derived from electric 
field correlations by considering the discontinuity of the normal electric field across the 
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surface of the plasma. The two-point electric potential correlations can be computed on 
a macroscopic scale by using linear response theory and macroscopic electrostatics. Let 
us put a test charge q at r. At some point r' the potential change 5$(r') due only to the 
plasma is related by linear response to the potential correlation at thermal equilibrium: 

5^{r') = -Pq < $(r)$(r') >^ , (2.1) 

where jS = {kBT)~^ with T the temperature and ks the Boltzmann constant, and 
< AB >^ means the truncated average < AB > — < A >< B >. 

The potential change 5$(r') can be computed using screening properties of the plasma 
and macroscopic electrostatic arguments provided that |r — r'| is large compared to the 
screening length. From now on we shall assume that this condition is satisfied. 

Let us apply this to the particular case of an insulated spheroidal plasma. We consider 
a ellipsoid of revolution around the z axis. Let 2b be its axial length and 2a its diameter. 
Let (P = 6^ — a^. Then \d\ is the distance between the foci of the spheroid. li cP > 
then the spheroid is prolate, otherwise it is oblate. We use spheroidal coordinates (^, rj, 0) 
defined by 

y^[{e-d^){l-v')]'^^sm<j>, (2.2) 
z = ^r]. 

The boundary between the plasma and the vacuum is then the spheroidal surface 
defined by 

e = 6. (2.3) 

Let us consider first the case where r and r' are inside the plasma. Due to the 
plasma's screening properties the charge g at r will be surrounded by a polarization cloud 
of microscopic dimensions carrying a charge —q giving to 5$(r') a contribution — g/|r — r'|. 
Since the plasma is insulated a charge +q spreads on the surface of the spheroid ^ = b. 



This surface charge gives to 5^{r') another contribution equal to q/C with the capacitance 
C given by [15] 



C = 



d 



(2.4) 



Qo(b/d) 

where Qo is a Legendre function of second kind. Using the total 5$(r') in (2.1) gives, for 
r and r' inside the spheroid 



P < $(r)$(rO >^= - i Qoib/d) . 

r — r d 



(2.5) 



If r' is outside but r is still inside the plasma the potential change at r' created by 
the surface charge q is equal to qQoii' / d) / d. So we have in that case 



/3 < $(r)$(r') >^ 



1 



r — r' 



d 



(2.6) 



If r and r' are both outside the plasma the potential change at r' is q[G{r, r') — 
|r — r'|~^] where G(r, r') is the potential at r' when a unit charge is put at r outside the 
insulated spheroidal conductor. G is given by [15] 



G(r,r') =-Y,i2n + l)Y,eU-iy 



n=0 



m=0 



(n — m)! 
(n + m)! 



-I 2 



X 



(2.7) 



+ 



Qo{^/d)Qo{^'/d) 



Qoib/d) 

where and are associated Legendre functions of the first and second kind respec- 
tively and €m = 2-Smo is the Neumann factor. The Pn{i/d)Q'^{^'/d) or P^ {i' / d)Q'^ {i / d) 
terms in (2.7) come from an expansion of |r — r'|~^, and finally the electric potential cor- 
relation is given by 



n 



(3 < $(r)$(rO >T=- ^(2n + 1) ^-(-1)' 



n=l 



m=0 



{n — my. 
{n + my. 



cos[m((/) — (f)')] 



X priv)prw)Qn{mQnie/d)^;}^^ 



(2.8) 



The surface charge correlation is 



1 



(2.9) 



where E^^'^°^^\r) denotes the hmit of the normal component of the electric field when r 
approaches the surface from the inside (outside). Using the expressions (2.5), (2.6) and 
(2.8) for the electric potential correlations we find 

1 P-(P d^G{r, r') 



P < (T(r)(j(r') > 

which finally gives 

P < a{r)a{r') >^= 



(2.10) 



r,r'£surface 



^ J] (2n + l)e. 



{ATr)H^{h^-dWW-d^rj'^) fi ^0 
in - m)\ , Q'^'ihld) 

X - — f7Pr(^)-pr(V) r^;,„.J cosm(0 - 00 . 



(2.11) 



(n + m)! " ^'^ " 

In the case of a globally neutral spheroid, < (j(r) >= and the truncation sign T may be 
omitted. 

From this last expression we can recover the charge correlation for some particular 
geometries. For example if d goes to zero, we have the case of a spherical plasma. In that 
limit h becomes the radius of the sphere and Q'^'{h/d)/Q'^{h/d) —{n + l)d/h then 
equation (2.11) becomes 



13 < (j{v)a{v') 



-I oo n / M 

1 V ^ V ^ , , , , (n — my. 

5^5^(2n+l)(n+l)6, 



Xir(^)^r(^')cosm(0-0')- 
The sum can be performed to give the already known result [2, 16] 



(2.12) 



(3 < a{r)a{r') 



87r2&3 



(2sinf)^ 



+ 



(2.13) 



where a is the angle between r and r'. 

Another special case is the cylindrical geometry obtained taking the limit b ^ oo, 
then a is the radius of the cylinder. In that case it is interesting to define k = n/b. The 
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sum over n times becomes an integral over k, b/d rsj I -\- rj rsj z/b and using the 
asymptotic expansions 

QniWd) ~ (inrKmika) , (2.14a) 

(2.14b) 



V nn 



in — m) kz 

^ ^2 



where the are modified Bessel function of the third kind, we find 

m=0 



1 °° k K' (ka) 

(3 < a(r)(7(r') >'^= -— V e„cosm(0-0') / - J^)^ [ cos k{z - z') dk . (2.15) 



III. COLLECTIVE MODES AND CORRELATIONS 

When the microscopic detail is disregarded, the thermal fluctuations are expected to 
be correctly described by the set of collective modes. If each collective mode n is considered 
as a harmonic oscillator of frequency a;„, the electric potential associated with this mode 
is of the form 

$n(r)e-*'^"* + '^J^e^'^'A /V2 , (3.1) 

with an amplitude of $n such that the corresponding average energy be /cbT, at temper- 
ature T. Then the time-displaced potential correlation will be 

< $(r, t)$(r', t') >^= Re ^ < '^Jr)^r.{r') > g-^'^^^"^') , (3.2) 

n 

and the correlation functions of the other functions can be deduced from (3.2) 

In the general case t ^ t' , expression (3.2) is expected to depend on the magnetic 
field applied to the plasma. However, the static limit t = should be magnetic field- 
independent. 

The explicit calculation of (3.2) for a spheroidal plasma in a magnetic field would 
involve complicated expressions which are not very illuminating. Therefore, only the spe- 
cial case (A) of a spherical plasma without a magnetic field is considered here. However, 
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as exercises, we consider simpler models, with magnetic field, on which it can be explic- 
itly checked that the static limit t = t' is field-independent. These models are (B), a 
plasma along a plane boundary in a magnetic field normal to the boundary, and (C), 
a two-dimensional plasma (with two-dimensional logarithmic Coulomb interactions) in a 
disk with a magnetic field normal to the plasma plane. 



We consider in this section a spherical one-component plasma of radius R composed 
of particles of mass m and charge g in a uniform charged background with charge density 
—qno, without a magnetic field. Experimentally this could be achieved in a Penning trap 
in the Brillouin regime (in the rotating frame the plasma behaves as an unmagnetized 
plasma) or in a Paul trap (as stated above, the confining fields play the role of the uniform 
neutralizing background) . 

We use spherical coordinates {r,6,(p). We compute the time-dependent correlations 
which can be seen as the sum of the contributions from the different collective modes each 
one oscillating at its own frequency. It should be noticed that the time-dependent corre- 
lations can also be computed by a generalization of the linear response method explained 
in section II, now using the dynamical linear response theory. 

The linearized equations of motion for the electric potential the volume charge 
density p and the current density j inside the plasma are 



A. Spherical plasma without magnetic field 



47rp, 



(3.3a) 



dt 
dp 

dt 




(3.3b) 



v.j. 



(3.3c) 



where ojr, = 



(47rq'^no/m)^/^ is the plasma frequency. We look for a mode of frequency uj. 
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Manipulating equations (3.3), we find for inside the plasma, the equation 

eA$ = , (3.4) 

where e = l — uj'^/uii'^. The equation for $ outside the plasma is the usual Laplace equation 

A$ = . (3.5) 

The problem has been reduced to the electrostatic problem of a dielectric filling the 
sphere of radius R. Equations (3.4) and (3.5) must be supplemented with the boundary 
conditions 



$ ^ when r +oo , (3.6a) 
lim $(r) = lim $(r) , (3.6b) 
edr^R', e, (j)) = dr^{R+, 9, 0) . (3.6c) 



Equations (3.4), (3.5) and (3.6) have two types of solutions: 

1) Surface modes: for e ^ 0, ^ satisfies the Laplace equation inside and outside the 
plasma. One finds modes (3.1) with 

^nmW - I ^^^2m^-/-lym(^^^)^ ii r > R, ^'^''^ 

and 

co^ ^ujI = — ivl , (3.8) 

" 2n+l P ^ ^ 

where n and m are integers (n > and |m| < n) and are the spherical harmonics. 

Equating the average potential energy of this mode to k-QT/2 gives the average squared 
amplitude 
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Using the time-displaced analog of (2.9), one find that the surface modes contribute to the 
time-dependent surface charge correlation < a{r, t)cr(r', t') >^ a term 



E 



(47r) 



2 < \An\^ > y-(^, 00 cosa;n(t - t') = 



(3.10) 



2) Volume modes: when e = 0, then a;^ = uj^. There is now an infinite number of 
modes for each (n,m): = and any ^'"^ = f{r)Y^{9, </>) with f{R) = is acceptable. 
However, a complete basis for the can be choosen as the eigenfunctions of the Laplacian 
with Dirichlet boundary conditions: 

$-(r) = 5]a,/,(r), (3.11) 

7 

with A/^ = A^/^, f^{R, 6', (^) = and / |/^(r)|2 dr = 1. 

Equating the average potential energy of each mode to A;bT/2, we find 
P < |a^|^ >= —Att/Xj. So, the contribution of the volume modes to < $(r, t)$(r', t') 
is 

J2 < 1^71' > fy{r)fy{r')cosiOp{t - t') = -AttUbT J] A; V^(r)/^(r') cos ujp{t - t') 

7 7 

= -47rA;BTGD(r, r') cos ojp{t - t') , 

(3.12) 

where Gd is the Green function of the Laplacian with Dirichlet boundary conditions on 
the sphere: 

1 r r ?? "1 

(3.13) 



r , i? 
— r r 

R r 



— r — r 



'1-1 



And finally the contribution of the volume modes to < a"(r, t)a{y\ t') is found to be 

kBT 1 



— 5- cosa;p(t — t') , 
87r2 (2i2sin(a/2))^ 

where a is the angle between r and r'. Putting (3.10) and (3.14) together 
(3 < a{v, t)a{v', t') ^ 1^^^^^' ^)WW:W) cos Un{t - t') 



(3.14) 



n,m 



+ 



(3.15) 



87r2 (2i?sin(Q;/2))^ 



cos a;p(t — t') . 



11 



For t = t' the sum in (3.15) can be performed and we recover the static result (2.13). 
B. Plasma in a half space with a magnetic field 

Let us consider now a one-component plasma filling the half space z < with a uniform 
magnetic field in the z direction, B = Bz. Sum rules for the time-dependent correlations 
have been obtained for this case using the dynamical linear response theory [17]. Here we 
use the collective mode method to find expressions for the correlations valid macroscopi- 
cally. This collective mode method has been used in the case B = in [12], here we extend 
it to the case B 0. With the same notation as in section III-A, the linearized equations 
of motion now are 

A$ = -47rp, (3.16a) 

^ = -/V$ + njAz, (3.16b) 
at 47r 

S = -V.j, (3.16c) 
ot 

where = qB /m is the cyclotron frequency. 

As above we look for a mode with frequency uj. From equations (3.16) we find for $ 
inside the plasma 

V . eV$ = , (3.17) 
where e is now the plasma dielectric tensor defined in Cartesian coordinates by 

(3.18) 




with ei = 1 — u;p/(a;^ — 1]^), €2 = l]a;p/[a;(a;^ — J]^)], and €3 = 1 — cu^/a;^. As to $ outside 
the plasma {z > 0), it obeys 

A$ = 0. (3.19) 
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We now have the problem of an anisotropic dielectric filling the z < half space. The 
boundary conditions are 



$ ^ when z +oo , (3.20a) 
^x,y,0-) = ^x,y,0+), (3.20b) 
esdM^, 2/, 0") = a,$(x, y, 0+) . (3.20c) 

Let us first consider the case ei/es < which give modes with a frequency in the ranges 
< < mm{u!p,Q) called magnetized plasma modes and max(a;p,f2) < |a;| < Qu = 
{ujp + n^)^/^ called upper hybrid modes. 
We look for a mode of the form 

^ (^e''"' + Be-'h^)e''^^-^^, for ^ < 0, . . 

"PkW - I + 5)e-/cx.e^k^-'-x, for ^ > 0, ^'^''^^^ 

where k = kj^ + k\\z, with kj^ in the xy plane, k± = |kj^| and rj^ = {x,y,0). Laplace 
equation (3.19) is satisfied for 2; > and equation (3.17) gives the dispersion relation 

1 - 41 ^11' + - ]k^' = 0. (3.22) 



Solving this equation we find two modes: one upper hydrid mode with frequency a;+ and 
one magnetized plasma mode with frequency a;_ : 



2 1 
^ 2 



(3.23) 

Equation (3.20c) gives a relation between the incident and reflected amplitudes A and 



B 



B 1 — 6163 

It should be noted that = 1: there is total reflection on the surface z = 0. 

The energy of this mode is a quadratic form of ^ and $. In the particular case B = 
this quadratic form is diagonal [12], which means that potential and kinetic energy have 
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equal averages and by the energy equipartition theorem each average is k^T jl (this was 
the case of section III-A). If S 0, the potential and kinetic average energies are different 
(take for example the simple case of a single charged particule with a circular trajectory in 
the plane normal to a magnetic field when the whole energy is kinetical). But we can still 
use the energy equipartition theorem and say that the total average energy of each mode 
is equal to k^T. 

The total average energy for a large volume V of plasma is found to be 

< E >= V — — ^ ^ . 3.25) 

Equating < > to /cb^ gives the average squared amplitude < \A\^ > of the mode. 

Finally these modes will give a contribution to the electric potential time-dependent 
correlation < $(r, t)^{r', t') >^ equal to 

V J [< $k,+ (r)$k,+ (r') > cosa;+(t - t')+ < $k,-(r)$k,-(rO > cosa;_(t - t') . 

(3.26) 

where $k,+ is the electric potential for the upper hybrid mode and 3>k,- the potential for 
the magnetized mode. 

To this contribution we must add the one from the possible modes in the range 
min(a;p, n) < |a;| < max(a;p,ri), called evanescent modes. In this case ei/es > and 
we look for a solution of the form 



(r) = I Ce^;;--+V^^-, if . < 0, (3.27) 

This form satisfies equations (3.17), (3.19) and (3.20a,b). Equation (3.20c) implies 63 < 
and €163 = 1. This means that we have an evanescent mode only if Q < u>p with frequency 
given by a;^ = a;^ = {fl'^+oj'^) /2. To compute < |C|^ > we proceed to compute the average 
total energy of the mode and equate it to k^T. In this case the energy is proportional to 
the surface S of the boundary between the plasma and the vacuum. We finally find 
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This mode adds a contribution to < $(r, t)$(r', t') >^ equal to 

< $kx(r)$kx(rO > cosa;e(i - ^0 • 



(3.29) 



In equation (3.26) it is convenient to make a change of variable in the integral over k\\ 
and have an integral over uj. In this way we can express the electric potential correlation 
< ^{r,t)^{r' ,t') >^ in terms of its Fourier transform C^^{k±, z, z' ,oj) with respect to 
time and the x and y coordinates: 

da;C**(kx,z,/,a;)e--(*-*')+^'^-('---'-'-). (3.30) 

C^^{'k±, z, z' ,u>) = if |a;| G]min(a;p, il), max(a;p, ri)[U]riu, +oo[ except maybe at 
CO = :hu>e if Q < u!p. And if co is not in that range 
-if 2; > and z' > 0: 

463 (-ei/e3)^/% _fe,(,+,.^ 



/?C$$(kj_, z, z',oj) 
-if 2; < and z' < 0: 

/3C<i>$(kj_, z, z' , iv) = 



u)k± 1 — €163 



(3.31a) 



a;A;_Le3(-ei/e3)V2 
1 



cos 



k^{-ei/esy/'{z-z') 



1 — €163 



(1 + 6163) cos \k±{-ei/e3)^/\z + z')\ (3.31b) 



+ 2e3(-ei/e3)'/' sin \k±{-ei/es)^/\z + z') 



-if 2; < and z' > 0: 



/?C$$(kj_, 2;, z'jLj) = 



1 e 



-k I z' 



ujk± 1 — 6163 



sm 



i-ei/e^y/'^k^z 



+ 63(-ei/e3)^/'cos [(-ei/e3)'/'/c±^] j 



(3.31c) 



In the case Q < ujp it must be added to /?C$$(kj^, z, z' ,u!) the term corresponding to 
the evanescent mode: 

exp —k±{z + z') if 2; > and z' > 0, 



(5(a;_a;,)+5(a; + a;J) J exp/c^^3f3^(^ + ^') if^<Oand^'<0, (332) 



exp /e±( ^2_Q2 -^ ~ -2') if -2 < and z' > 0. 
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From these expressions we can compute the surface charge correlation, using equation 
(2.9). The surface charge correlation Fourier transform is found to be 

4A;±(1 - €3)2 

(3.33) 



'^e3(— ei/e3)V2(i _ ^163) 



+ 



47r ojI - 



The last term is to be included only when Q < u>p. 

Now, let us briefly show how these results can be obtained using dynamical linear 
response. If we put at 2; = 2;' an oscillating charge density 5{z — z') exp z(kx • rj^ — ojt), the 
electric potential change at z is x(kj^, z, z\lo) expi(kj^ • rj^ — Loi). The response function 
X is related to the Fourier transform C$$(kj^, 2;, 2;', a;) of the time-dependent correlation 
< $(r,t)$(r',f) in the non-perturbed system by the fluctuation-dissipation theorem: 



^ , 1 , 
/3C<|,$(k_L, 2, 2 ,0;) = Imx(k_L, 2,2 ,0;), 

TTO; 



(3.34) 



and x(k_L, 2, 2', uS) = \E'(k_L, 2, 2', u) — 2n exp[—k±\z — z'\]/k± where \E' is the total electric 
potential, due to the plasma and the external charge, solution of 



V . eV 



*(k_L, z, z', u))e^^^-^'''-^^ = —4:7r6{z — z') exp z(k_L • rj_) , 



(3.35) 



where e is the dielectric tensor given by (3.18) if 2; < or equal to 1 if 2; > 0, and the 
boundary conditions (3.20). 

Solving equation (3.35) gives for x 
-if 2; < and 2;' < 



X(k^,2;, z'.uj) 



27r 
'kl 



1 /i^ 



£3 V ei 



+ 



exp 



ei 



1 + esV^i/^s 
exp[— /cj_|2 — z\\ 



exp 



-/CiJ-(2 + 2') 
63 



(3.36a) 
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-if 2; > and z' < 



/ X 27r 
k± 



exp 



1 + esV^i/es 
— exp[— — z'W 



-k±z + ^1 — k±z' 



(3.36b) 



-if z > and z' > 



X0i±,z,z',uj) = ^-^^^^^^exp[-k±{z + z')] . (3.36c) 

From these expressions it is easy to verify that equation (3.34) leads to (3.31) and 
(3.32) and therefore both methods give the same results. Furthermore, for t = t', equation 
(3.34) and Kramers- Kronig relation 

"+°°ImxO^±,z,z',io) 



7rRex{^±,z,z',0)=P f 

J — c 



(3.37) 



give the well-known static correlation [2] 
/3 < $(r)$(r') >'^= 



r — r 



if 2; < or ^' < 0, 



[\r±-r^'\^ + {z + z')^] if 2 > and^' > 0, 



which is, as expected, independent of the presence of the magnetic field. 



C. Plasma in a disk with a magnetic field 



(3.38) 



Here we consider the model of a two-dimensional a one-component plasma in a disk 
of radius R with a magnetic field normal to the plane where the disk lies. The particles 
interact through the two-dimensional Coulomb potential — Inr. We look for modes with 
frequency oj. The dielectric formalism is also valid in this case. Only the definition of the 
plasma frequency is slightly changed: ujp — {2nq'^no/my^'^. The equation for the potential 
$ inside the disk is 

eiA$ = 0. (3.39) 
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while outside the disk it is the usual Laplace equation. In polar coordinates (r, 9) the 
boundary conditions are 



$ —5- when r — > +oo , (3.40a) 
lim $(r) = lim $(r) , (3.40b) 

ei i€2 — = . (3.40c) 

or roB or 

As in section III-A, there are two types of solutions: 

1) If ei 7^ 0, ^» satisfies the Laplace equation inside and outside the disk 

if^<^' f3 41) 

^ ' ~ 1 A^i^2|-lr-l-le^™^ if r > i?, > 

and there are two possibles frequencies for each integer m (m ^ 0), u) = sgn[m)u!±, where 
sgn(m) denotes de sign of m and 

uj± = (^-n ± + 2ujfj /2 , (3.42) 

Equating the total average energy of the mode to k^T gives 

/5 < \Am? >= I , (3.43) 

where u is the frequency of the mode and u' the other root of (3.42). The contribution 
from these modes to < a{r,t)a{r' ,t') >^ is 
\ \ im(9-e') 

/cbTV— ^-^ ^^g-isgn(m)a;_(t-t') _^_e-»sgn(m)a;+(t-t')| (3.44) 

^ (7ri?)2(a;+ -a;_) L J 

2) If ei = 0, a;^ = a;^ + ri^ and any $ satisfies equation (3.39) inside the disk. However, 
writting $ as a Fourier series in e*"^^, boundary conditions (3.40) and Laplace equation 
for $ outside the disk implies that $ = outside the disk. The situation is similar to the 
one in section III-A 2). $ can be written in the base of the eigenfunctions of the Laplacian 
with Dirichlet conditions on the boundary. Then equating the average total energy of each 
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mode to ksT and following the calculations from section III- A 2) we find the contribution 
from these modes to < $(r, t)$(r', t') 



-27r/cBr— i-GD(r, r') cos + Q^/^t - t')] , (3.45) 

where Gd now is the Green function of the Laplacian with Dirichlet boundary conditions 
on the disk: 

GD(r, r') = — [In |r - r'| - In |r - {R/r'fr'l] . (3.46) 
27r 

In two dimensions there is another mode corresponding to a; = 0. The dielectric formalism 
does not applies here (e2 diverges when a; = 0). Dealing directly with the equations 
of motion, we can show that $ = outside the disk. Then the contribution of this 
mode can be computed in the same way as in the case ei = 0. The contribution to 
< $(r, t)$(r', t') >^ is found to be 

-2nkBT— -GD(r, r') . (3.47) 



Putting together all contributions gives 



P<a{r,t)a{T',t') >^ = J] 



me 



im(9-9') 



(7ri?)2(a;_l_ — U- 



^^+^1 8 (Tripsin f)' ■ 

(3.48) 

For t = t' , the sum in equation (3.48) can be performed and the static result [2], indepen- 
dent of the magnetic field B is recovered 

(3 < a{r)a{r') ^ . (3.49) 

8 (Tripsin |) 

CONCLUSION 



Since one-component plasmas have been obtained experimentally, it would be inter- 
esting if the static or dynamical correlations could be measured and compared to the 
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expressions obtained here. It should be noted that our calculations give only the dom- 
inant term of the expansion of the correlations in powers of the distance r between the 
points. It has been shown that there are other terms in the asymptotic expansion [18] 
of the time-dependent correlations (behaving like for the potential correlations in an 
infinite plasma). These algebraic corrections vanish in the static case. 
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